Quasiuniversal connectedness percolation of polydisperse rod systems 
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The connectedness percolation threshold (n c ) and critical coordination number (Z c ) of systems 
of penetrable spherocylinders characterized by a length polydispersity are studied by way of Monte 
Carlo simulations for several aspect ratio distributions. We find that (i) rj c is a nearly universal 
function of the weight-averaged aspect ratio, with an approximate inverse dependence that extends 
to aspect ratios that are well below the slender rod limit and (ii) that percolation of impenetrable 
spherocylinders displays a similar quasiuniversal behavior. For systems with a sufficiently high 
degree of polydispersity, we find that Z c can become smaller than unity, in analogy with observations 
reported for generalized and complex networks. 

PACS numbers: 64.60. ah, 61.46.Fg, 82.70.Dd 



Idealized elongated objects such as perfectly rigid 
cylinders, spherocylinders and prolate spheroids are pro- 
totypical models for a wide array of technologically rel- 
evant systems that include liquid crystals, nanocompos- 
ites based on filamentous fillers as well as fiber-reinforced 
materials. Percolation phenomena involving dramatic 
increases in, e.g., structural rigidity and electrical and 
thermal conductivities of composites with increasing filler 
loading are currently of particular interest [l|. These in- 
creases are caused by the formation of an infinite cluster 
of in some sense connected particles at the critical load- 
ing, i.e., the percolation threshold. 

It has been established by analytical [IHII and numeri- 
cal P4TlT| studies that for dispersions of sufficiently elon- 
gated objects of identical size and shape, i.e, "monodis- 
perse" objects, the geometric percolation threshold ex- 
pressed in terms of the critical volume fraction of parti- 
cles is inversely proportional to the aspect ratio of the 
filler particles. This property is exploited in the fab- 
rication of conducting polymeric composites with very 
low conducting filler contents. Depending on the pro- 
duction processes of the composites, however, the filler 
particles almost invariably exhibit a pronounced poly- 
dispersity in both size and shape [lj, [H|. Although 
it represents a possible factor behind huge quantitative 
discrepancies between theory and experiments [lil [Tfij . 
such polydispersity has received relatively little atten- 
tion in terms of theoretical modeling until fairly recently 
[E EH S3 j Achievement of a theoretical understanding of 
how the continuum percolation of fibrous fillers is affected 
by polydispersity is thus key to the controlled design of 
a large class of composite materials for practical particle 
size and shape distributions. 

Recent analytical results obtained from integral equa- 



tion methods [161 ] and from an heuristic mapping onto a 
generalized Bethe lattice [TtJ predict that in the slender 
rod limit, where the particles have asymptotically large 
values of the aspect ratio, the volume fraction at the per- 
colation threshold is inversely proportional to the weight 
average L w = (L 2 )/(L) of the rod lengths, where the 
brackets imply number averages over the distribution of 
rod lengths L. This Letter presents Monte Carlo (MC) 
results for the percolation threshold of isotropically ori- 
ented spherocylindrical particles with length polydisper- 
sity and having aspect ratios ranging from ~ 1 to several 
hundreds. We show that the percolation threshold of 
polydisperse, interpenetrable spherocylinders is a nearly 
universal function of L w over the entire range of aspect 
ratios considered. In addition, the percolation thresh- 
old closely follows the predicted 1/L W behavior even for 
particles with aspect ratios that are considerably smaller 
than the slender rod limit, thus generalizing the current 
theory. 

For systems of impenetrable spherocylinders with fixed 
\J (L 2 )/D, where D is the diameter of the hard core, we 
show that the percolation threshold is nearly indepen- 
dent of the length distribution. Finally, we find that the 
critical coordination number per particle at the percola- 
tion threshold (denoted Z c ) can be smaller than unity 
for polydisperse systems. Although similar observations 
have been reported for a number of complex networks and 
in systems of hypersphcrcs in high-dimensional spaces 
[l8l ] , this finding is novel in the context of the continuum 
percolation of three-dimensional objects. 

We generate isotropically oriented distributions of pen- 
etrable rods by randomly placing N penetrable shero- 
cylindcrs with a distribution of lengths L and identical 
diameter S within a cubic box with periodic boundary 
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FIG. 1: (Color online) (a) critical distance ratio R between 
polydisperse and monodisperse spherocylinders as a function 
of the fractional occupancy p of the longer rods with shorter 
rod length fixed at Li = 20 in units of 5 c o = 2/irp{L 2 ) (see 
text). The solid lines represent Eq. (J2j) . (b) The critical 
distance ratio R in units of Ro calculated from Eq. © for 
L1/L2 — 3 as a function of the length Li of the shorter rods 
and for selected values of p. 



conditions and side length C. As a measure of the concen- 
tration of the spherocylinders we shall use the dimension- 
less density r) = p(v), where (v) = (tt/6)S 3 + (tt/A)S 2 (L) 
is the number-averaged volume of the spherocylinders, 
(L) = JdLLf(L) is the mean (number-averaged) rod 
length for a given distribution f(L) of lengths, and 
p = N/ C 3 is the number density of the particles [19| . 

We consider two spherocylinders as being connected 
if they overlap geometrically. The percolation threshold 
is identified by ascertaining the minimal diameter S c (for 
fixed value of p) for which a cluster of connected particles 
spans the entire cubic box. This definition is equivalent 
to the usual procedure of finding a critical density p c 
of spherocylinders with fixed diameter, and has the ad- 
ditional advantages of: (i) being computationally more 
convenient, and (ii) allowing a more direct relation to 
the conductivity a of rods through the critical distance 
approximation a cx exp(— 2<5 c /£), where £ is the tunneling 
decay length ■ 

In the following we shall use the critical distance 
<5 c o, defined as S c q = 2/7rp(L} 2 , as our unit of length 
for polydisperse systems of penetrable spherocylinders. 
This quantity corresponds to the critical distance ob- 
tained from the second virial approximation formula 
rj c = (1/2)S c q/L for the critical concentration of a system 
of monodisperse spherocylinders with identical lengths 
L 3> <5 c o chosen to coincide with (L). 

To find 5 C we employ the clustering method described 
in Ref. [2fjl which allows computation of the spanning 
probability as a function of the sphero cylinder diameter 
S for fixed density p (Supplemental Material [21]). Fig- 
ure [TJa) shows the results obtained for polydisperse sys- 
tems with a bimodal length distribution f(L) = pS(L — 
Li) + (1 - p)S(L - L 2 ) with L 2 = 20 and L x > L 2 , 



where < p < 1 is the number fraction of long rods. In 
the figure we display the ratio R (symbols) of the crit- 
ical distances for the polydisperse rod system to those 
for monodisperse systems of spherocylinders with lengths 
equal to (L) = J dLLf(L), which for the particular dis- 
tribution considered corresponds to (L) = pLi + (l—p)Li. 
The ratio R of the critical distances is systematically re- 
duced by polydispcrsity and displays a minimum that 
becomes deeper and moves towards smaller values of p 
as L1/L2 is increased, implying that a small fraction 
of longer rods can substantially lower the percolation 
threshold. 

This trend is in full agreement with the theory of 
Ref. [TH based on the second virial approximation to the 
connectedness Ornstcin-Zcrnike equation, which predicts 
for (L)/S c > 1: 



2 L lt 



(i) 



where L w ~ (L 2 )/(L) is the weight-averaged rod lengths. 
From this equation, and by using r] c ~ p(tt/A)5 2 (L) for 
{L)/8 C ^> 1, the critical distance is predicted to follow 
S c = 2/ttp(L 2 ). The reduction factor Ro — S c /5 c q pre- 
dicted by the theory is thus 



Ro — 



(L) 2 = (pLjLi + l-p) 2 
(L 2 ) p(Li/L 2 )2 + l-p' 



(2) 



where the second equality applies for the bimodal length 
distribution. As shown in Fig. [lja) the MC findings for 
R are in semi-quantitative agreement with Eq. ([2]) (solid 
lines), although R is consistently slightly smaller than Rq. 
This discrepancy could arise from the circumstance that 
the values of L\ and Li used in the simulations may be 
insufficiently large to achieve the slender rod limit, which 
is a prerequisite for the validity of Eq. ^ . We examine 
this issue in Fig. [TJb), which shows R/Rq as a function 
of L2 (up to Li = 150 in units of (5 c o) for L\jLi = 3 and 
selected values of p. For Li > 20, our MC results are 
less than 10% smaller than R/Rq = 1. Furthermore, for 
Li > 50, R/Rq appears to increase monotonically (albeit 
somewhat slowly), which may indicate that the slender 
rod length limit R/Rq = 1 could ultimately be reached 
(for any p) only for very long rod lengths. 

The MC results shown in Fig.[T]and the relatively small 
deviations from Eq. © suggest that, for rods with identi- 
cal radii, L w is the key quantity that controls the percola- 
tion threshold for mixtures of penetrable spherocylinders. 
This is demonstrated in Fig. [5] where r\ c is shown as a 
function of L w /S c for various bi-disperse (open symbols) 
and monodisperse (+ signs) systems of spherocylinders 
(in the latter case L w is identical to the unique particle 
length) . Results for systems of spherocylinders for which 
the lengths follow Weibull and uniform distributions are 
shown in Fig. [5] by filled circles and squares, respectively 
[2TI HU . Surprisingly, all of our data collapse onto a sin- 
gle curve over the entire range of L w /S c > 1, implying 
that rj c is a quasiuniversal function of L w /S c independent 
of the particular distribution considered. 
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FIG. 2: (Color online) Critical reduced density r\ c as a func- 
tion of L W /5 C for monodisperse (plus symbols), bi-disperse 
(open symbols), Weibull (filled circles), and uniform (filled 
squares) distributions of spherocylinder lengths [2l], [22] . The 
solid line represents Eq. {TJ. Inset: p<5j? as a function of 
\J (L 2 ) /5 C calculated for mododisperse and bidisperse impen- 
etrable spherocylinders with hard core diameter D. The x 
symbols are the results for both the monodisperse and poly- 
disperse penetrable rods of the main panel. 



This finding is rather unexpected because the observed 
quasiuniversality extends well below the slender rod limit 
of Eq. ([I} (solid line), which is approached by the MC 
data to within less than 10% only for L W /5 C > 200. Fur- 
thermore, even though Eq. ([1]) might be expected to ap- 
ply only asymptotically for L W /5 C 3> 1, we observe that 
the data for L W /8 C > 10 are well fitted by a(L w /S c )~° 
with a = 0.165 ± 0.009 and /3 ~ 1.080 ± 0.002. The in- 
verse scaling of rj c with L w thus applies approximatively 
even for spherocylinders with very modest aspect ratios. 

We have also examined the effects of length polydisper- 
sity on the percolation of impenetrable spherocylinders 
with identical hard-core diameters D. Two impenetrable 
spherocylinders {D ^ 0) arc considered to be connected 
if their surfaces approach closer than 5. The percola- 
tion threshold for a given density p of the particles is 
identified, as before, by the critical distance S c . In the 
slender rod limit the percolation threshold of impenetra- 
ble rods is predicted to follow (j> c = D 2 /{25 C L W ) [H,[I3, 
where cf> c ~ p(ir/4)D 2 (L) is the critical volume fraction 
for the hard-core particles. By noting that Eq. (QJ is 
the percolation threshold for penetrable rods, and since 
■q c ~ p(tt/4:)S 2 (L}, we see that for sufficiently elongated 
rods the percolation relation 

P 5l = (2MS 2 /(L 2 ) (3) 

is predicted to be satisfied by both hard and penetrable 



rods, independent of their length distribution. 

We have generated by MC simulations equilibrium dis- 
persions of impenetrable spherocylinders with different 
length (L) distributions. The inset of Fig. [2] shows p5% 
as a function of \J (L 2 ) /6 C for monodisperse systems with 
L/D = 10 and 20 (filled symbols) and for two bidisperse 
cases with L\, L 2 , and p chosen as to give \J (L 2 ) /D = 10 
and 20 (open symbols) [2l| . Although for computational 
reasons the rod lengths considered by us are not large 
enough for our results to fulfill Eq. ([3]), we see neverthe- 
less that for a given yj (L 2 )/D, pS^ is essentially indepen- 
dent of the particular rod length distribution. Further- 
more, the calculated p8\ values for increasing J (L 2 ) /D 
tend to follow the same functional behavior of the inter- 
penetrable spherocylinders (x signs in the inset of Fig. [5]) 
[231 1. This latter feature suggests that for sufficiently 
large yj (L 2 )/D there exists a universal relation of the 
form p5c = F{yJ (L 2 ) /S c ), which is expected to reduce to 
Eq. ([3]) for yj (L 2 ) /5 C 3> 1 and that applies to both pen- 
etrable and interpenetrable spherocylinders over a wide 
range of yj (L 2 )/S c values. 

Although currently there is no theoretical explanation 
for the quasiuniversal dependence reported in Fig. [2j 
a partial understanding may be achieved by follow- 
ing the method developed in Ref. [Hj]. In this for- 
malism, applied here for simplicity to penetrable rods, 
the overall cluster size S satisfies S = (T(L))l where 
T(L) - p{C+(L,L',5 c )T{L')) L , = 1, and C+(L,L',6 C ) 
is the orientation-averaged connectedness direct correla- 
tion function at zero wave vector. Within the plausible 
ansatz C+(L,L',S C ) = LL'c u + (L + L')S 2 c w + S^c 00 [3], 
where the coefficients {Cy} are assumed to depend only 
upon the packing fraction, it is found that for systems 
with different length distributions but equal values of 
L w j 'S c , S diverges at percolation thresholds that differ by 
~ a 2 5 c /L w for L w /S c > 1 and by - (L w /S c )a 2 /(I + a 2 ) 
for L W /5 C < 1, where a 2 = (L 2 )/(L) 2 - 1 is the scaled 
variance. Since the scaled variances for all length distri- 
butions considered in this work were always smaller than 
~ 50% [IJ, for L w /S c > 10 the expected deviation from 
universal behavior is thus only ~ a 2 5 c /L w < 5%, which 
is consistent with the results of Fig. [5] Although we ex- 
pect that systems with values of a 2 much larger than 
those considered by us would imply a stronger deviation 
from universality, a 2 < 0.5 is nevertheless representative 
of the scaled variances observed in several real polydis- 
perse systems of rod- like particles [HI, 13 • 

The quasiuniversal dependence of the percolation 
threshold upon L w implies a general non- universality of 
the critical coordination number Z c , where Z c denotes 
the average number of contacts per rod at the percola- 
tion threshold. This is best viewed for the case of ran- 
domly placed and oriented overlapping objects for which 
Z c = »7c(^ex)/ where (u C x) is the excluded volume av- 
eraged over the orientations and the rod lengths. Given 
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FIG. 3: (Color online) Critical coordination number Z c as a 
function of L W /8 C for polydisperse and monodisperse sphero- 
cylinders. The symbols have the same meaning as in the main 
panel of Fig. [2] 



that rj c depends on L w /S c , while 

M , {{L)/5 c f 



l + (3/2)(L)/r5 c 



(4) 



depends on the rod lengths through (L) /S c , we see that 
mixtures of rods with equal r\ c (i.e., equal L w ) may have 
rather different Z c if the distribution f(L) of the rod 
lengths is such that (L) ^ L w . 

Figure [3] shows the critical average number of connec- 
tions per rod calculated from Z c — r] c (v cx )/(v) for the 
same mixtures of penetrable spherocylinders considered 
in Fig. [2] We have verified that the configurational av- 
erage of the connection per rods at S c coincides with the 
excluded volume formula, as expected. We see that in 
general, Z c is sensitive to the extent of polydispersity, al- 
though in the limit L w /S c — > it can be expected that Z c 
should coincide with the result for identical, overlapping 
spheres, namely Z c ~ 2.74 [25J. In particular, Z c is al- 
ways larger than unity and approaches Z c — > 1 asymptot- 
ically in the slender rod limit for monodisperse systems. 
In contrast, for distributions with sufficiently large val- 
ues of the variance and of L w / 5 C , polydisperse systems of 
rods may display fewer than one connection per particle 
at the threshold, i.e., Z c < 1. 

This latter feature is somewhat novel since the contin- 
uum percolation of objects randomly dispersed in a three 
dimensional space is usually characterized by the condi- 
tion Z c > 1 [T^]. Indeed, to the best of our knowledge, 
percolation occurring with Z c < 1 has been reported only 
for penetrable identic al hy perspheres in spaces of dimen- 
sionality exceeding 12 [18[ and in random or complex net- 
works that are not embedded in a physical space. For ex- 
ample, given an uncorrelated network with nodes having 



a distribution of coordination numbers z, upon random 
removal of nodes the network becomes disconnected at a 
critical node occupation probability p c = (z)/((z 2 ) — (z)) 
26 , [27j , which results from the irrelevance of closed loops 
28[. The critical coordination number Z c = p c (z) is thus 



Z r .= 



{zf 



(5) 



which can be smaller than unity when the node degree 
distribution is such that (z 2 )/(z) — (z) > 1. 

The results of Fig. [3] that show that polydisperse rod 
mixtures may display Z c < 1 suggest that these systems 
may relate to such classes of generalized graphs that can 
exhibit the same feature [2{|. Indeed, as shown in Ref . Ii~7l. 
Eq. (J5]) is also the critical coordination number of a gen- 
eralized Bethc lattice that by construction lacks closed 
loops. Hence, by applying the mapping (z) — > 2(L)/S C , 
(z 2 ) — > 4(L 2 )/5 2 formulated for polydisperse slender rods 
in Ref. ujl we find 



(Vf 
(L 2 ) 



< 1, 



(6) 



which is qualitatively consistent with the behavior of Z c 
seen in Fig. [3J A physical explanation for the observa- 
tion that Z c < 1 for sufficiently polydisperse rod systems 
is provided by the fact that the percolating cluster is 
predominantly comprised of the longer rods in the sys- 
tem, and the shorter rods have a greater likelihood of 
being isolated [17| . An interesting corollary arising from 
this interpretation is that, as in generalized random net- 
works where targeted removal of highly connected nodes 
enhances the percolation threshold (2(1 [27| , preferential 
removal of the longer rods from the system may lead to 
similar enhancement of the critical concentration. 

In conclusion, we have studied by MC simulations the 
effects of length polydispersity on the percolation thresh- 
old of penetrable spherocylinders. We find a quasiuni- 
versal dependence of the percolation threshold on L w 
that extends well below the slender rod limit considered 
in Refs. 0, [l?}- For systems of impenetrable sphcro- 
c ylinde rs we find that universality is fulfilled for a given 
\J (L 2 )/D, where D is the hard-core diameter. The pre- 
dicted quasiuniversality could be tested by experiments 
in systems of conducting fibrous fillers by altering the 
distribution of the rod lengths, e.g., by sonication, and 
measuring the resulting change in the percolation thresh- 
old. Furthermore we have demonstrated that the average 
number of connections per rod at the percolation thresh- 
old can be smaller than unity for random distributions of 
rods that are sufficiently slender and polydisperse. This 
finding reveals an intriguing analogy with the case of ran- 
dom percolation in complex networks. 

B. N. acknowledges support by the Swiss National Sci- 
ence Foundation (Grant No. 200020-135491). 
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SUPPLEMENTAL MATERIAL 
A. Calculation of the critical distance 

For both penetrable and impenetrable spherocylinders 
we follow same the method to calculate the critical dis- 
tance S c . Namely, for a given number density p of sphe- 
rocylinders which are either penetrable or impenetrable 
with hard-core diameter D, we coat each spherocylindcr 
with a penetrable shell of thickness 6/2, and we consider 
two spherocylinders to be connected if their penetrable 
shells overlap. For penetrable systems (i.e., for D = 0) 6 
represents the diameter of the penetrable spherocylinder. 
For each realization of the system, we compute through 
the clustering method described in Ref. [l[ the minimum 
value of (5 such that a cluster of connected spherocylin- 
ders spans the entire sample. By counting the number 
of instances that sample-spanning clusters appear for a 
given <5, we construct the percolation probability curve 
P(S). 

As in the main text, for systems of penetrable rods 
we adopt as unit of length the quantity S c o = 2/7r ( o(L} 2 , 
which corresponds to the critical distance in the second 
virial approximation for monodispersc spherocylinders 
with length fixed at (L). Examples of P(S) obtained 
from 500 realizations of polydisperse (solid lines) and 
monodisperse (dashed lines) systems of penetrable sphe- 
rocylinders are shown in Fig.|4]for different numbers N of 
spherocylinders with density fixed at p = 7.89 x 10 -4 . For 
the polydisperse cases we have considered a bi-disperse 
length distribution f(L) = pS(L - L x ) + (1 -p)5{L - L 2 ) 
with L\ = 60, L2 = 20, and the number fraction of long 
rods p = 0.21. The monodisperse systems were generated 
by spherocylinders of length equal to (L) = J dLLf(L), 
which for the particular distribution considered corre- 
sponds to (L) = 28.4. Figure |4] reveals that the spanning 



1 

poly mono 




- N 


= 2.5 10 3 l/ t 


- N 


= 4 10 3 / / 


N 


= 1 1 o 3 / / 


_ N 


= 20 10 3 // // 




// il 






0.02 


0.03 0.04 0.05 




8/<L> 


Percolation 


probability P(S) for different numb 



N of penetrable spherocylinders as a function of 5/{L) for 
monodisperse (dashed lines) and bi-disperse (solid lines) sphe- 
rocylinders. The parameters of the bi-disperse distribution 
are Li = 60, L2 — 20, and p = 0.21. The monodis- 
perse systems are generated by considering rods with iden- 
tical lengths coinciding with (L). All lengths are in units of 
Sco = 2/np(L) 2 . In the figure the number density is fixed at 
p = 7.89 x 10" 4 for all cases. 



probabilities for the bi-dispersed systems are shifted to 
lower values of <5 when compared to the P(S) curves for 
the monodisperse case, indicating that the polydisperse 
systems percolate at smaller volume fractions. 

For both the polydisperse and the monodisperse cases 
the curves for the two highest values of N intersect at 
approximately P = 1/2, which we take as our criterion 
for identifying the critical distance 5 C . For the particular 
case of Fig. H we find S C /{L) ~ 0.046 and 0.032 for the 
monodisperse and polydisperse cases, respectively. 

Our results for penetrable spherocylinders shown in 
Figs. 1, 2, and 3 of the main text have been obtained by 
considering simulation box sizes C such that C/L\ > 5, 
where Li is the largest rod length for any given distri- 
bution, and the number N of particles exceeds 2 x 10 4 . 
The resulting critical distances have been obtained by 
adopting the criterion P{8 C ) = 1/2. 

Figure[5]shows the spanning probability P(S) obtained 
from 300 equilibrium configurations of bi-disperse sys- 
tems of impenetrable spherocylinders with L 1 /D = 30, 
L 2 /D = 10 and p = 3/8 and for different values of the 
number density p. From the largest to the lowest den- 
sities the number N of spherocylinders decreases from 
N = 7000 to N = 3000 and the box size C increases from 
C ~ 4Xi to C ~ lOLi. We have used a fitting to a simple 
sigmoidal function to evaluate the critical distance from 
P(dc) = 1/2 and from the mean of the distribution func- 
tion dP(5)/d8. The two methods give values of 5 C which 
differ at most by a few percent. 



B. Bi-disperse, Weibull, and uniform distributions 

To avoid exceedingly large computational times, we 
have been careful to choose rod length distribution func- 
tions with length Li of the longest rod not exceeding 
about 4 times that (L2) of the shortest one. Despite 
of this constraint, we still have been able to generate 
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FIG. 5: Percolation probability P(5) for different densities p 
of impenetrable and bi-disperse spherocylinders as a function 
of S/D, where D is the hard-core diameter. The parame- 
ters of the bi-disperse distribution considered in the figure 
are Li/D = 30, L2/D = 10, and p = 3/8. The corresponding 
value of sJ(L 2 )/D is 20. 
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FIG. 6: Scaled variance a 2 = (L 2 }/(L) 2 — 1 as a function of the number fraction p of long rods for the bi-disperse distribution 
of lengths used in the calculations. 
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FIG. 7: Discretized Weibull (left panel) and uniform (right panel) distribution functions of the spherocylinder lengths. 



length distributions with large scaled variances a 2 — 
(L 2 )/(L) 2 — 1. Among the different distributions con- 
sidered, the bi-dispcrse one, i.e., f(L) = pS(L — Li) + 
(1 — p)8(L — L 2 ) with < p < 1, had the largest a 2 for 
a given L\jL 2 [2[. Figure [6] shows a 2 = p(l — p)(n — 
l) 2 /[p{n - 1) + l] 2 where n = Li/L 2 = 2, 2.5, 3, and 4 as 
a function of p. We see that for p ~ 0.2 and L1/L2 > 3, 
a 2 is well above 30%. 

In Fig. [7] we show the discretized Weibull (left panel) 
and uniform (right panel) distribution functions of 
the rod lengths used in our study on the penetrable 
polydisperse spherocylinders. The discretized Weibull 
distribution is defined as f{L{) = exp[— (Li/\) k ] — 
exp[— (L i+ i/X)/ k ], where Li = i (i = 1,2,3,...) arc 



the rod lengths We have used k = 6 and A = 3, 
15, 60, and 110. The corresponding scaled variance 
(j 2 = (L 2 )/(L) 2 — 1 is g 2 — 4%. To guarantee that the 
ratio of lengths between the longest and shortest rods 
never exceeded ~ 4, the distribution was truncated (and 
subsequently normalized) by eliminating from the sam- 
pling all spherocylinders with f{Lj) smaller than 10~ 2 of 
the maximum of the distribution. 

Uniform distributions of rods (right panel of Fig. [7]) 
have been constructed from f(L) = l/(£i — L 2 ) for L 2 < 
L < L\ and f(L) = otherwise, with L\jL 2 = 4 and 
L 2 = 1, 5, 20, and 50. For all cases the scaled variance 
is er 2 = 12%. 
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